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Introduction: Chemically and microstructurally 
complex altered rims around grains in the finest size 
fraction (<20 (J.m) of the lunar regolith are the result of 
multi-stage processes involving both solar ion radiation 
damage and nanoscale deposition of impact or sputter- 
derived vapors [1], The formation of the rims is an 
important part of the space weathering process, and is 
closely linked to key changes in optical reflectance and 
other bulk properties of the lunar surface [2]. Recent 
application of field-emission scanning transmission 
electron microscope techniques, including energy dis- 
persive X-ray spectral imaging, is making it easier to 
unravel the “nano-stratigraphy” of grain rims, and to 
delineate the portions of rims that represent Radiation- 
Amorphized (RA) host grain from overlying amor- 
phous material that represents vapor/sputter deposits 
[2]. For the portion of rims formed by host grain amor- 
phization (henceforth called RA rims), we have been 
investigating the feasibility of using Monte Carlo-type 
ion-atom collision models, combined with experimen- 
tal ion irradiation data, to derive predictive numerical 
models linking the width of RA rims to the grain’s in- 
tegrated solar ion radiation exposure time. 

Model Approach: Several in-situ TEM ion irradi- 
tion studies have precisely calibrated the critical ion 
doses associated with ion radiation-induced amorphiza- 
tion in rock-forming silicates [3,4,5]. For total ion en- 
ergies corresponding to interaction between the ions 
and target atoms by predominantly nuclear/elastic col- 
lisions, silicates typically show abrupt amorphization 
when the total nuclear/elastic collision energy depos- 
ited into a unit volume of crystal reaches a critical 
level, irrespective of the type or energy of the incident 
ion [3,4], The critical deposited energy E d for amorphi- 
zation can serve as a normalizing parameter that allows 
experimental radiation-induced amorphization data 
obtained at one set of ion total energies to be used to 
predict the ion doses for amorphization at other ion 
energies. The computational tool that allows the two 
energy regimes to be compared is provided by a Monte 
Carlo-type ion-atom binary collision code such as 
SRIM [6]. Because SRIM also provides data for depos- 
ited collision energy as a function of depth, it ulti- 
mately allows the depth of the amorphous layer in an 
irradiated target to be calculated as a function of in- 
creasing ion dose, as calculated from grain exposure 
time and time-averaged solar wind ion flux. 


Experimental Inputs: We determined the critical 
ion doses and E d for amorphization for a pure synthetic 
enstatite (MgSiOV) and San Carlos olivine (F090) using 
the in-situ IVEM-Tandem TEM ion irradiation facility 
at Argonne National Laboratory. The details of the 
irradiation procedures are described in [3,4,7]. For 
irradiation with 1 MeV Kr ions, an energy regime 
where nuclear/elastic collisions dominate (similar to 
lower energy solar wind H + and He ++ ), we determined 
E d (pyroxene) = 3400 eV/nnr 3 and E d (olivine) = 3000 
eV/nm 3 . The olivine value is in excellent agreement 
with previous in-situ irradiation studies [3]. The pyrox- 
ene value has to our knowledge never been previously 
determined by in-situ irradiation means. 

Model Inputs: The experimental E d values define 
the critical level of deposited ion collision energy at 
which pyroxene and olivine become amorphous, inde- 
pent of the irradiating ion and its energy, as long as 
ion-atom interactions are in the nuclear stopping re- 
gime. Because this requirement holds for ions at typical 
solar wind energies, our experimental E d values are 
applicable to modeling solar wind radiation damage. 
Using the SRIM Monte-Carlo code we calculated the 
combined total nuclear collision energy deposited by 
solar wind H + and He ++ as a function of depth in py- 
roxene and olivine for various grain exposure times in 
the top-most layer of the lunar soil. The E d values were 
used to determine the maximum depth of amorphiza- 
tion for a given exposure time from the corresponding 
total ion dose and deposited energy versus depth curve. 
In this approach, one of the main variables is the time- 
averaged flux versus energy distribution for the solar 
wind ions. For the current results we used solar wind 
speed variation data from [8], combined with time- 
averaged solar wind flux estimates that we have previ- 
ously made for lunar grains [9]. Other variables in- 
clude the time-average distribution of ion incidence 
angles, and the presence or absence of sputtering. For 
the former we ran SRIM with a plasma-like 271 stera- 
dian distribution of ion incidence angles. For the latter 
we simulated sputtering by using SRIM’s ability to 
incrementally change the starting depth of the incident 
ion. 

Model Results: For equivalent solar wind irradia- 
tion conditions, olivine and pyroxene have nearly 
equivalent curves for the growth of their radiation- 
induced amorphous rims as a function of exposure age 




Figure 1: Modeled widths of solar ion-amorphized 
rims on lunar soil grains as a function of exposure time 
in the top-most layer of the lunar soil. Results are for 
solar wind speeds narrowly distributed about 450 
km/s. 

(Fig. 1). The “plateau” feature of the curves is typical 
and reflects the drop off in deposited energy versus 
depth as ions reach the end of their range combined 
with the decrease in solar wind flux at solar wind 
speeds above the 850 km/sec range (e.g., 3.7 keV H + 
energy). Pyroxene shows slightly higher (~5 nm) pla- 
teau values for amorphous rim widths relative to oli- 
vine for long exposure times (Fig. 1). This slight dif- 
ference is unlikely to be detectable in the population of 
lunar grains given the effect of more important factors 
such as variable grain exposure times. When sputtering 
is added to the simulation using lunar sputtering rates 
suggested by Kerridge [10], 24 nm of surface removal 
over 100,000 years results in a moderate 9% decrease 
in the RA rim width (Fig. 1). Overall, sputtering has 
the net effect of decreasing the total nuclear/elastic 
energy deposited at a given depth relative to the un- 
sputtered case, extending the exposure times needed to 
attain the same depth of amorphization. However, our 
initial results indicate that sputtering rates much higher 
than those suggested by Kerridge [10] would be needed 
to produce detectable differences in RA rim widths 
across a given population of grains. The Figure 1 
curves show the growth of amorphous rims for a very 
narrow solar wind flux versus energy distribution cen- 
tered around a mean solar wind speed of 450 km/s. 
Curves for a more realistic set of solar wind speed dis- 
tributions [8], based on calculations for pyroxene are 
shown in Figure 2. The plateau values for the rim 
widths for typical grain exposure times of 5000- 
150,000 years agree well with typical rim widths meas- 
ured by TEM [1,2]. This result confirms the reason- 
ableness of the model assumptions and approach, an 



Figure 2: Modeled widths of solar ion-amorphized rims 
on pyroxene grains in the lunar soil as a function of expo- 
suree time for different solar wind speed distributions. 


important first-level goal of the current study. Also 
shown are changes in the curves that result from mod- 
erately skewing the solar ion flux-versus-energy distri- 
bution to both lower and higher energies relative to a 
normal (Gaussian-like) distribution (Fig. 2). Such 
changes are part of year-to-year and longer timescale 
variations in solar activity [8]. The growth curves show 
little change for a negatively skewed distribution 
largely because of the form of the He ++ stopping power 
curve at lower energies. But this same factor makes the 
model output sensitive to a positive skew, suggesting 
that RA rim widths may help identify lunar soil grains 
exposed during extended periods of higher solar activ- 
ity, assuming the exposure ages of these grains can be 
independently estimated. 
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